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ABSTRACT 

We present a new scheme for rotations of a charge qubit associated with a singly 
ionized pair of donor atoms in a semiconductor host. The logical states of such a qubit 
proposed recently by Hollenberg et al. are defined by the lowest two energy states of 
the remaining valence electron localized around one or another donor. We show that an 
electron located initially at one donor site can be transferred to another donor site via 
an auxiliary molecular level formed upon the hybridization of the excited states of two 
donors. The electron transfer is driven by a single resonant microwave pulse in the case 
that the energies of the lowest donor states coincide or two resonant pulses in the case 
that they differ from each other. Depending on the pulse parameters, various one-qubit 
operations, including the phase gate, the NOT gate, and the Hadamard gate, can be 
realized in short times. Decoherence of an electron due to the interaction with acoustic 
phonons is analyzed and shown to be weak enough for coherent qubit manipulation being 
possible, at least in the proof-of-principle experiments on one-qubit devices. 
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1. INTRODUCTION 

Solid-state systems are of great interest in the search for a scalable quantum computer 
technology. Several schemes for solid-state quantum information processing have been 
proposed [HIS El- For example, the coherent control of superconducting qubits |1] and 
their coupling |5| have been demonstrated, the qubits being encoded in the states of 
a Cooper-pair box. One promising area of current investigation is concerned with the 
semiconductor-based devices. In the Kane proposal JBJ, the qubits are defined by long- 
lived nuclear spins of phosphorous dopants in a silicon host. They are manipulated by 
external surface gates and RF magnetic fields. While long coherence times of nuclear spins 
make the Kane scheme very promising, the single-spin measurement remains a significant 
challenge [7]. This also concerns an alternative Si:P architecture that uses electron spin 
states as qubits jS]. 

Along with spin-based qubits, the charged-based qubits in semiconductors are currently 
discussed as well. The logical states of a charge semiconductor qubit may be formed 
by, e. g., the ground and excited states of the electron in the single quantum dot jT] 
or the spatially separated states of the electron in two different quantum dots jHl US 
rn\ IT^ ITB] . In spite of the fact that decoherence of the charge-based qubits is rather 
strong imiini, the charge qubits are nevertheless believed to be realizable at the present 
technological level due to their short operation times JHl- One of the obstacles to the 
practical realization of scalable quantum computation in the system of quantum dots is 
that it is extremely difficult, if at all, to manufacture a set of quantum dots with identical 
or at least predetermined characteristics each. This complicates the issue, introducing the 
errors into the operations with qubits ^7] and generating a need for numerous ancillary 
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corrective gates. In this respect, it would be more reasonable to make use of natural 
atoms (instead of "artificial" ones) as the localization centers for the electrons carrying 
the quantum information. Recent advances in manipulation with single atoms on the 
solid surface PH] and atomically precise placement of single dopants in semiconductors 
[ini I2ni naake possible the construction of rather complex solid-state atomic architectures. 

Recently, HoUenberg et al. proposed a two-atom charge-qubit scheme JHI and reported 
the first results on its fabrication and characterization 20J for the case of phosphorous 
dopants in silicon. In that scheme, the buried donor charge qubit consists of two dopant 
atoms ~ 50 nm apart in a semiconductor host. One of the donors is singly ionized. The 
logical states are formed by the lowest two energy states of the remaining valence electron 
localized at the left or the right donor, |0) = \L) and |1) = \R), see Fig. 1. The qubit is 
controlled by the surface electrodes through adiabatic variations of the donor potentials. 
Initialization and readout of the qubit are facilitated by a single electron transistor. The 
coupling of such qubits via the Coulomb interaction allows, in principle, to realize the 
conditional two-qubit gates [TBj . 

It was shown in Ref. |T^ that although the coherence time Tcoh ~ 1 ns for charge-based 
qubits is much shorter than for their spin-based counterparts, the corresponding gate 
operations times are also shorter, of order Top ~ 50 ps. Note, however, that the ratio of 
Top/ Tcoh ~ 10^^ seems to be insufficiently small for the fault-tolerant scalable quantum 
computation being possible [21]. Here we propose an alternative scheme for operations 
with buried donor charge qubits, instead of applying biases to the surface gates. Our 
scheme is based on the effect of electron transfer between the lowest states localized at 
different donors upon the influence of a resonant pulse ^ or two resonant pulses [22] • 



3 



Such a transfer occurs via an excited molecular level of the double-donor system and 
allows for implementation of different one-qubit rotations. The operation times can be 
made orders of magnitude shorter than in the original proposal |1.6| . 

The paper is organized as follows. In Section II, we describe a three-level model for the 
resonant electron transfer between the donors and briefly discuss the relevant one-electron 
states of a P2' molecular ion in Si. Next we present the analytical solution for the unitary 
electron evolution under the influence of microwave pulses. In Section III, we show that in 
the P^:Si system it is possible to realize various one-qubit operations, including the NOT 
gate, the phase gate, and the Hadamard transformation. Decoherence due to the electron 
interaction with acoustic phonons is studied in Section IV. Discussion of the results is 
given in Section V. 



We consider a singly ionized pair of phosphorous atoms embedded in silicon. The 
remaining valence electron is described by the Hamiltonian 



where En and \xn) are, respectively, the one-electron eigenenergies and eigenstates of the 
molecular ion Pg'iSi. In general, to calculate the energy spectrum and the wave func- 
tions {r\xn) of the single-electron/double-donor system beneath the surface, one should 
account for the conduction-band anisotropy, the inter- valley terms, the surface effects, the 
potentials induced in the substrate by the gate voltages, etc. This necessarily requires 
numerical calculations, see, e. g. Ref. [22] ■ We note that although the conduction-band 
edge of bulk Si has six degenerate minima, it has been shown both experimentally [23 



II. MODEL FOR THE RESONANT ELECTRON TRANSFER 



-^0 = X! En\Xn){Xn 




n 
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and theoretically that substitutional impurities break the translational symmetry of 
the crystal lattice, thus lifting the degeneracy. The spacing between the energy levels in 
the ground-state and excited-states multiplets may be further increased through appro- 
priate choice of the gate potentials. Anyway, to quantify the structure of the P2^:Si energy 
spectrum and wave functions, one should make sophisticated numerical calculations for a 
specific donor configuration. In this paper, however, we restrict ourselves to a semiquan- 
titative consideration based on an isotropic effective mass approximation [211 that allows 
for an explicit analytical solution. Then the problem reduces to that for the hydrogen-like 
molecular ion with the effective Bohr radius ^ 3 nm and the effective Hartree unit of 
energy E* = e'^/ea*^ ~ 40 meV, where e = 11.7 is the dielectric constant for silicon fI7\ . 
The energy spectrum of the H.2 ion for different atomic separations is known with high 
accuracy [2H1- 

We approximate the Hamiltonian Hq in Eq. (QJ by the reduced three-level Hamiltonian 

Hr = Ei\xi){xi\ + E2\X2){X2\ + Etr\xtr){xtr\ , (2) 

where and \x2) are the lowest molecular states, Isag and 2pau, whose wave functions 
are, respectively, symmetric and antisymmetric about the midpoint of the line joining the 
two donors, see Fig. 2, and \xtr) is one of the excited molecular states discussed below. 
It is convenient to go from the states and IX2) delocalized over the P2^:Si ion to the 
states \L) = [\xi) + IX2)]/ and |-R) = — |x2)]/\/2 localized at the left and the right 
donor, respectively. For donor separations R» a*^, the wave functions (r|L) and {r\R) 
are almost indistinguishable from the one-electron Is-orbitals of the corresponding donor 
atoms. 

The states \L) and \R) form the qubit logical states |0) and |1), respectively. These 
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states are well defined if the thermal energy ksT is much lower than the differences 
AE^i = E^ — Ei and AE32 = E3 — E2 between the energy E^ of the excited molecular state 
Ixs) and the energies Ei and E2, respectively. At i? >> one has Ei ^ E2 ^ —E*/2 
and E3 ^ -E*/8, so that A^si ^ AE32 ^ 3^78 ^ 15 meV. Since the states \L) and \R) 
are not the exact eigenstates of the Hamiltonian Hr, in the absence of external fields the 
initial qubit state |\&(0)) = a;|L) + will evolve with time as 

= e-'^^*/''|^(0)) = e-'^^*/^||^(0))+z(/5-a)e-'^^2^*/2^sin(AE2iV2;i) \L)-\R) 

(3) 

where AE21 = E2 — Ei. Note that at t << to = %/ AE21 the initial qubit state remains 
almost unchanged (not counting the common phase). Since dX x = R/ a*^ » 1, the value 
of AE21 is exponentially small [23 1^0] , 



AE. 



21 _ A^„-X-l 



Axe '-^ 



(4) 



E* 

the period to ~ ti/AE2i it takes for the qubit state to change is rather long, > 1 fis at 
i? > 60 nm. In what follows we shall consider the processes taking place in time intervals 
much shorter than to and hence ignore the off-diagonal term [— (Ai?2i/2)|L)(i?| + h.c] in 
Hr that gives rise to the electron tunneling \L) ^ \R). Then the Hamiltonian (0) takes 
the form 



y El + E2 



\L){L\ + \R){R\ 



Etr\xtr){xtr\ , (5) 
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where {Ei + E2)/2 ^ Ei ^ E2 at R » a*^. In the general case that the qubit is biased 
by the gate voltages, the energies El and Eji of the lowest states localized, respectively, 
at the left and the right donor differ from each other. In this case, the localized states 
are all the more good approximations to the energy eigenstates, and the Hamiltonian Hr 



reads 



Hr ^ EL\L){L\+En\R){R\ + ETR\xTR){XTB^ 



(6) 



Now let the buried donor charge qubit interact with an external electromagnetic field 
E(t). Then the Hamiltonian becomes 



with = {xtr\ — er\L) and = {xtr\ — er\R) being the electric dipole moments for 
the transitions \L,R) ^ \xtr) between, respectively, the localized states \L) and \R) and 
one of the excited molecular states \xtr) delocalized over the double-donor system. For 
definiteness, we choose this state to be the third one-electron state Ixs) of the molecular 
ion P2':Si. At = Er and R/a*^ > 6, this is the 3dag state whose wave function {r\x3) 
is symmetric about the midpoint of the line joining the two donors and has its maxima at 
the donor locations [21], see Fig. 2. If the donors are arranged along the x-axis, the state 
Ixs) is formed upon the hybridization of \2S)l,r and \2Px)l,r atomic states of the donors, 
and (r|x3) in the vicinity of the left/right donor equals to [{r\2S)L,R T {^I'^Px) l,r]/2 at 
R » a%. Note that for such a choice of the state Ixs), the electric field should have a 
nonzero x-component in order that /j ^ 0. 

We consider two cases: (a) Ei = Er ^ Ei and (b) E^ ^ Er, the desired value of the 
difference Er — E^ being discussed below. In the case (a), we suppose E(t) to oscillate 
at a frequency u = {Etr — El^r)/?i, 



H{t) = Hr + V{t) 



(7) 



where the interaction term V{t) is 



V{t) = E{t) dL\XTR){L\ + dR\xTR){R\ + h.c. 




E(t) = Eo(t) cos(cjt) 



(9) 
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where Eo(t) is the slowly varying envelope of the field. Making use of the resonant 
approximation i. e., omitting the rapidly oscillating terms with the frequencies 

±{uj + ExR/h — Ei^ji/h) from the Hamiltonian, we have 

1 



V{t) = -e-^-* 



XLit)\xTR){L\+XRit)\xTB){R\ 



+ h.c., (10) 



where Xifiif) = Eo(t)dL^/j. In the case (b), the field E(t) has two components oscillating 
at frequencies ujl = {Etr — EL)/h and ur = {Etr — ER)/h, 

E(t) = Eoi(t) cosKt) + Eo2(t) cos{ujRt + 0) , (11) 

where is the phase shift between the two components. In the resonant approximation 
[32] one has 

V{t) = le-'^^'XL{t)\xTR){L\ + \e-'^'''-''f'XR{t)\xTR){R\ + h.c. , (12) 

where XL,R(t) = Eoi,2(i)dL,_R- In this paper, we restrict ourselves to the rectangular pulse 
shape, so that Eo(t) in Eq. and both Eoi(t) and Eo2(i) in Eq. ifTT]) are constant at 
< t < Top and zero elsewhere. 

It is straightforward to solve the non-stationary Schrodinger equation for the state 
vector |^(t)), 

^^^^^ = H{t)\m) , (13) 

with the Hamiltonian H(t) in Eq. ((Tj) given by Eqs. (0) and (fTUI) in the case (a) or ® and 
(fT^ in the case (b), and to find the coefficients Ciit), CR^t), and CxRif) in the expansion 
of \^{t)) in terms of the states and \xtr), 

\<^{t)) = Ci(t)e-^^^*/"|L) + C«(t)e-^^-*/"|i?) + C™(t)e-^^™*/^|xTi?) , (14) 



provided that |^I/(0)) = a\L) + I3\R), where |ap + = 1. In the case (a) we have 



CL{t) = a 
Cnit) = - 



2\Xl 


to 


|Al| 




Ar 


2 



2ArA 



a- 



Al|2 + |a^|2 

2|A^|2 



|Al|2 + |A^|2 



sin^(nt) , 



n i4.\ ■ + ^-^^ • ^on+^ 



v/|AlP + |Ar| 



(15) 



where 



AlP + IA^I 



Ah 



(16) 



In the case (b), the coefficients Ciit), Cji{t), and CxRit) are also given by Eqs. (fT3j) 
and (fTBj) with the only exception that \r should be replaced by Xrc''"^. From Eqs. (fT3|) 
and ()16|) one can see that at t = Top = 7Tk/2Q (hereafter is a positive integer) the 
coefficient Ctr vanishes, so that the state vector |^(t)) remains in the qubit subspace 



{\L),\R)} and |Cz.(rJ|2 + |C«(r, 



op J 



1. In particular, if Ci(0) = 1 and Ci{(0) = 0, 



then Cl{top) = and Cr{top) = ±1 at A^^ = =fA_r and odd k, i. e., there is a complete 
population transfer \L) \R), see Ref. 0. So, the auxiliary excited state \xtr) plays the 
role of the "transport" state, in that it assists the qubit evolution by means of the electron 
transfer between the states \L) and \R) as the pulse is on but remains unpopulated after 
the pulse is offi 



III. QUBIT ROTATIONS 

In this Section we show that the auxiliary-state-assisted electron transfer between the 
two donors allows for various qubit rotations. In the case (a) that the two donors in the 
molecular ion P2 :Si are equivalent, i. e.. El = Er and |Al| = \Xr\, the qubit state |^'(t)) 



9 



at the operation time Top remains unchanged, 



(17) 



if Top = nk/Q, or changes into 



P\L) + a\R) 



if Top = 7i{2k — 1)/2Q and = T^r, see Eqs. (fT^ and (fT3j) . The latter corresponds to 
the quantum NOT operation. 

The case (b) seems to be more reahstic because of the different local atomic surroundings 
of the donors in the pair due to both the uncontrollable damage of the host upon ion 
implantation and the probabilistic variations in the path taken through the substrate by 
each implanted ion ^Uj. Besides, the surface gates can be used to intentionally tune Ei 
and Eji to the predetermined values. Moreover, one can change the values of and Xr 
separately through the changes in the electric field amphtudes Eqi and -Eo2- It follows from 
Eqs. (fT^ and (fT3j) that the relative phase shift operation is implemented at Top = nk/Q, 



a\L) + (3e 



-i(ER-EL)Top/h 



\R) 



(19) 



while the value of Top = TT(2k — l)/2f2 corresponds to realization of the quantum NOT 



operation. 



P\L) + a\R) 



(20) 



if Al = tXr and = Trn+ {ER — EL)Top/2h (hereafter n is an integer), or to the Hadamard 
transformation. 



V2 



V2 



(21) 
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if {Eji — E^jTop/Ti = 27rm (where m is a positive integer). Here the plus sign corresponds 
to the values of = 27m and A^^ = — A/j(v^ — 1) or = 71(211 + 1) and \l = ^r{V^ — 1), 
and the minus sign corresponds to the values of (p = 27in and = ^r{V2 + 1) or 
(j) = n{2n + 1) and Al = -Xr{V2 + 1). 

So, various one-qubit operations can be implemented on the buried donor charge qubit 
through appropriate choices of the pulse frequency, phase, amplitude, and duration. Let 
us estimate the value of the operation time Top ~ 1/Q ~ ^/IAl^/jI ~ h/ea*^EQ, see Section 
II. For the field amplitude i?o ~ 1 V/cm one has Top ~ 1 ns. Increase in the pulse intensity 
will cause the value of Top to decrease down to the picosecond time scale, so that the value 
of be made orders of magnitude shorter than the period to it takes for the qubit 

state to change due to the direct electron tunneling \L) ^ \R), see Section II, as well as 
the operation times in the case that the qubit is manipulated by adiabatically varying 
the potentials of the surface gates [inj. Note that in the case (b) the energies El and 
Eji should be sufficiently different from each other in order all these operations could be 
implemented in short times to avoid decoherence, as discussed below. For example, at 
Top ~ 1 ps one should have Eji — El 3 meV. 

IV. DECOHERENCE EFFECTS 

The uncontrolled interaction of the quantum system with its environment leads to 
entanglement between the states of the system and the environmental degrees of freedom. 
This disturbs the unitary evolution of the system and results in the loss of coherence. 
There are various sources of decoherence in solids. For the charge qubit considered in this 
paper, the decoherence due to the phonon emission/absorption processes was studied in 
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Refs. jini I2E] and found to be much weaker than the decoherence due to both Nyquist- 
Johnson voltage fluctuations in the surface electrodes and 1/f noise from the background 
charge fluctuations. Note, however, that there are two mechanisms of the phonon-induced 
decoherence which are caused by either the energy relaxation processes or the virtual- 
phonon dephasing processes. Which one of those mechanisms is dominant, depends on 
the specific parameters of the quantum system and its environment, as well as on the 
operation times. Here we show that the dephasing processes play a decisive role in limiting 
the fault tolerance of the buried donor charge qubit. For simplicity, we consider the qubit 
at zero temperature and assume isotropic acoustic phonons with the linear dispersion law, 
cUq = sq, where s is the speed of sound. 

First we recall some general concepts concerning the transition probability for an elec- 
tron moving in the time-dependent potential. If an electron, being initially in the state 
\i) of the discrete energy spectrum, interacts with the harmonic field 



then the probability amplitude to find it in the state |/) at a time t is given by the 
following expression that results from the first-order perturbation theory 



where cUj/ = {Ei — Ef)/h. The common approach is to ignore the first term in Eq. (j2SI) 
and make use of the expression 





(23) 




(24) 
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thus arriving at the so called Fermi golden rule for the transition probability, 

h [uJif — ujy ri 

(25) 

where Ti^f{ijj) is the time-independent transition rate. The 5-function reflects the energy 
conservation, TiuOif = hu, for such a transition. 

Electron-phonon coupling in confined systems is described by the Hamiltonian 

Hei-ph = T.K^)pi^)[K + K] ' (26) 
q 

where 6q and are, respectively, the operators of creation and annihilation of a phonon 
with the wave vector q, /)(q) = / dre^'^'^ p{r) is the Fourier transform of the electron 
density operator p(r) = X^mn ^m(^)^n('^)l"^)('^l; '^(^) the microscopic electron- 
phonon interaction matrix element that can be expressed in terms of the deformation 
potential D and the density of the crystal p as 



with fl being the normalizing volume. If the harmonic field ()22p is associated with a defor- 
mation phonon having the frequency cijq, then, taking into account that the deformation 
fields produced by the phonons with different wave vectors are not correlated, one has for 
the total transition rate 



where 

F,f{c^) = X{ci){i\e^^^\f) . (29) 
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A. Decoherence during adiabatic variations of the surface gate potentials 

In the case that the buried donor charge qubit is controlled by the surface gates [15] , so 
that the state vector |\E'(t)) remains in the qubit subspace {\L), \R)} during the operation, 
and the overlap {L\R) is negligibly small, the Hamiltonian (j^U]) can be written in the spin- 
boson form |Sn], 

Hel-ph = &z Y.9i^) [bq + Lq] , (30) 
q 

where a, = \L){L\ - \R){R\ and 

A(q) 



{Lle'^^'lL) - {R\e'''''\R) . (31) 



2 

Since (r|L, R) = (vr(a|j)^)~^/^ exp(— |r — r^^ijl/a^) for Is-orbitals, where Tl^r are the donor 



coordinates, one has 



sm{q,R/2) 

[1 + {qa%r/4y 



9{^) = -a(q)^^^-— — 2 , (32) 



where Qx is the component of the phonon wave vector along the line joining the two donors, 
and we chose the origin of the coordinates in between the donors, so that = ^{R/2)ex- 
Fedichkin and Fedorov ^3] have shown that at T = decoherence upon implementing 
the phase operation emerges as pure dephasing, the electron density matrix being given 
by the general expression jSZl EH] , 



Pll(O) p^^{Q)e-B\^)+^iE,-E,)t/n 
PijL(0)e-^^W-(^«-^-)*/'^ p^,^(0) 



(33) 



with the spectral function 



There is no relaxation in this case since in order the phase operation could be implemented, 
the energies El and should be sufficiently different from each other J3], so that the 
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basis {\L), \R)} coincides with the energy basis of the electron in the double donor system, 
and the electron term (jU)) commutes with the interaction term (|3U|) in the Hamiltonian. 
As a result, the diagonal elements of the density matrix remain unchanged. On the other 
hand, decoherence upon implementing the quantum NOT operation (where El = En and 
the energy basis of the electron is formed by the states |xi,2) = \\L) ± \R)]/ v^, see Section 
II) was suggested to be caused by relaxation fT^, so that both off-diagonal and diagonal 
elements of the density matrix decrease exponentially with time, the relaxation rate r2^i 
being [11112^1, see Eq. (1^ . 

(35) 



where g2i = AE2i/sh, see Eq. (jH). 

Note, however, that the approximation ()25|) for Wi^f{uj,t) and, accordingly, the equa- 
tion (j^Hj) for Ti^f are valid provided the time t is sufficiently long, see Eq. (j211). To 
quantify the applicability of this approximation, let us analyze the more general expres- 
sion for Wi^f{t) that follows from Eq. (j2Sl), 

sin^ ( ^^M^^t] 

w^wW = 72El^^/(q)r . ^ I. J ■ (36) 

One can roughly distinguish two phonon contributions to Wi^f{t), one being from the 
"resonant component", i. e., from the (5-f unction-like peak of sin^ C^'-''^ '^'^ t) / {uif — cjq)^ as 
a function of q a.t q = qif = Uif / s, with the height and the width ~ 1/st, and another 
from the remaining " non- resonant background" of the phonon spectrum. The former can 
be estimated as 

W^fit) ^ ^\F,f{q.f)\H , (37) 
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and the latter as 

wBfit) - ^\F.f{qma.)\' (38) 

at Qif « Qraax and 

^S/W~^fMV./(w)P (39) 

at qif » qmax, where Qmax is the wave vector at which the function |Fjj(g)p has a 
maximum, and Aq is a characteristic width of |Fjj(g)p in the maximum. The specific 
values of Aq, qmax, and Fif^q^ax) depend on the specific type of wave functions (r|z) and 
(r|/) in the matrix element {i\e^'^^\f). Now if, e. g., qif « qmax and we are interested 
in the transition probability Wi^f(t) at a moment of time t such that sqfj\Fif{qif)\H « 
Aq\Fif{qmax)\'^, then wj;^j{t) << wj-^f{t), and hence the Fermi golden rule appears to 
be broken l4Uj . This is due to violation of the energy conservation at short times j^S]- 
An inspection of the phonon-induced transitions between the states \xi,2) = [\L) i: 
\R)]/^/2 of the double donor system with El = Er and the donor separation R » a*^ 
(these transitions are relevant for decoherence during the implementation of the NOT 
operation ^3]) provides an illustrative example of the departure from the Fermi golden 
rule. In this case (2|e^'i'"|l) = [(L|e^^''|L) - (i?|e^^'"|i?)]/2, so that F2i(q) = ^(q), see Eq. 
(1221), and the resonant component of the transition probability is W2^i{t) ~ q^^D'^t / phs^ , 
in accordance with the value of the relaxation rate r2-»i given by Eq. ()35p. Since the 
value of q2i = AE2i/hs decreases exponentially with R, see Eq. ()H), the value of r2^i 
decreases exponentially as well, going below 10^ s~^ at R/a*^ > 10, see Fig. 5 in Ref. 
PB] . On the other hand, since g2i << Qmax ~ V^b, we have W^2^i('^) ~ D"^ / phs^ {a*^Y 
from Eq. More accurate calculations result in W^^X^it) = B^{t)/2, see Eq. (jSlD- If 

the operation time Top is long compared to the phonon transit time, a*^/s (~ 0.3 ps for 
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P^:Si), one has from Eq. (El), see Ref. |Tlj, 

so that the spectral function (jH^ appears to be a material constant, being about 6 ■ 10^'^ 
for the phosphorous donors in silicon ^3], where D = 3.3 eV, s = 9 • 10^ cm/s, and 
p = 2.33 g/cml Hence, W^^^iit) » W^Si(i) at R/a*^ = 10 and t«i^3 - 10"^ s, the 
time t being exponentially longer for larger values of R/a*^, and in any case longer than 
the operation time Top, see Section III. 

So, contrary to suggestions [11112^1 that phonon-induced decoherence in the case of the 
NOT operation is determined by the value of the relaxation rate r2^i given by Eq. 
we see that at sufficiently short operation times, decoherence in the cases of both phase 
and NOT operations is determined by the same spectral function -B^(t), see Eq. (j34j) . 
The distinction between the two cases is that the diagonal elements of the density matrix 
remain unchanged in the case of the phase operation since there is no relaxation, while 
they decay exponentially (along with the off-diagonal matrix elements) in the case of the 
NOT operation [H. 

B. Decoherence during the auxiliary-state-assisted operations 

Since the excited "transport" level \TR) becomes temporarily populated during the 
resonant-pulse operations on the P^:Si qubit, the phonon-induced electron transitions 
\TR) ^ \L,R) and \TR) ^ \xi,2) can have a detrimental effect on the qubit evolution, 
along with the transitions \L) ^ \R) and ^ \x2) studied above. Let us clarify what 
type of the phonon-induced electron transitions ("resonant" or " non- resonant" ) is domi- 
nant in this case and estimate the transition probability. We follow the line of reasoning 
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outlined above and start with calculations of the matrix elements {TR\e'^'^^\L,R) . For 
our choices of the "transport" state \TR) = Ixs) and the double donor orientation, see 
Section II, at R » one has \TR) ^ [\2S)l - \2P^)l + \'2S)r + |2P^)r]/2, where 
{r\2S)L,R = (87r(a|j)3)-V2(i _ |r _ ri_«|/2a|j) exp(-|r - rL,R\/2a%) and {r\2P,)L,R = 
(327r(a^)^)^^/^(x — xl,b) exp(— |r — rL^R\/2a*^). Neglecting the exponentially small over- 
lap between the localized atomic-like orbitals centered at different donors, we have 

{TR\e'^'^\L, R) = 2v/2^^^^^^id(^e^^«-^/2 . (41) 

[I + 

Since, depending on the relative values of and Er, the lowest energy eigenstates of 
P2 :Si are either \L) and |-R) (if E^ 7^ Er) or xi and X2 (if = Er), in order to find the 
probability WxRit) of electron escape from the "transport" state at T = 0, one should 
add up the probabilities of, respectively, \TR) \L) and \TR) \R) or \TR) \xi) 
and \TR) — > IX2) electron transitions. If the value of Er — E^ is much less than the 
difference between Etr and El^r, then in both cases we have the same result, so that 
WxRit) is given by Eq. (jHBj) . where now Uij ^ AE^i/h ^ 3E*/8h does not depend on R 
at R >> a*^, and 

|F.Kq)r = 16^^f^^±l^^|A(q)r (42) 

Taking into account that g^/a^ = Uija^/s ~ 3e^/8ehs ~ 8 >> qmaxCi*B ~ 1; it is 
straightforward to derive from Eq. ()36|) the following expressions for the probabilities of 
the "resonant" and " non- resonant" transitions, respectively, 

^^^^^^ ^ vrp;..2(«*^)3(*/«B) t:— —t (43) 

and 

176^2 



^^^^^^ ?,U^T^^phs\a%f{qifa%Y ' ^^^^ 
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It follows from Eqs. ^ and (gH) that W^/^it) = TrRt, where Ttr ^ 3 ■ 10^ and 
W^^{t) ^ 10^^, so that the "resonant" transitions are dominant at t > 0.3 ps. 

Now let us check what states out of those involved in the auxiliary-state-assisted qubit 
evolution are most sensitive to phonon-induced decoherence. As we have seen above, 
decoherence of the low-energy states \L) and |_R) (or and IX2)) is quantified by the 
error rate ^1], i. e., the error generated during the operation time, D{t) = B'^{t)/2 
3 ■ 10^^. This value is greater than W^)i{t) but less than Wj'^(t) at t > 10^"'^° s, where the 
processes of the spontaneous phonon emission by an electron temporarily occupying the 
"transport" level become prevailing. So, at Top < 100 ps, the error rate does not exceed 
the value of D{Top) ^ 3 ■ lO^^. 

V. DISCUSSION 

Fast auxiliary-state-assisted evolution of the double-donor charge qubit driven by the 
resonant electromagnetic field allows for implementation of various one-qubit rotations 
in very short operation times Top < 100 ps, thus minimizing the unwanted decoherence 
effects. At such times, the error rate due to acoustic phonons is D{Top) ~ 3 ■ 10^'^ at T = 0. 
At finite temperatures, such that fc^T > Tiujq, where huo = hs/a*^ ^ 2 meV for dephasing 
processes and "non-resonant" emission/absorption transitions, and huo = \Ei — Ef \ for 
the "resonant" \i) ^ |/) transitions, the error rate increases by a factor of ~ ksT /hujo. 

The most strong increase in the error rate at T 7^ occurs if the two donors in the 
molecular ion P2^:Si are equivalent since in this case the energies and Efi of the lowest 
localized states \L) and |-R) are equal to each other, and the difference E2 — E1 between the 
eigenenergies of the two lowest delocalized molecular states and IX2) is exponentially 
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small at large donor separations, e. g., E2 ~ Ei ^ 10~^ meV at i? = 60 nm, see Eq. 
(jH). To weaken decoherence, it would be reasonable to make use of the surface gates in 
order to increase the difference Eji — El up to E^ — i?^ ~ 1 meV so that the energy basis 
of the electron be formed by the states \L) and |-R) instead of the states and \x2)- 
In this case, the electromagnetic field should have two components driving the electron 
transitions |L) ^ \TR) and |-R) ^ \TR) between the states \L,R) and the auxiliary 
"transport" state \TR). 

At T 7^ 0, the processes of the phonon absorption by an electron temporarily occupying 
the "transport" state also contribute to decoherence. For our choice of the "transport" 
state, \TR) = Ixs), the state nearest to it in energy is the state Ix^)- In the case that the 
two donors are equivalent and R/a*^ > 15, this is the 4/cr„ state Ixi) ~ [|25')i — |2P^)i — 
1 25)/? — \2Px)r]/2 whose wave function (r|x4) is antisymmetric about the midpoint of 
the line joining the two donors At x = R/a*^ » 1 the energy separation [^01 

E4^ ~ E3 = E*{x^/4:) exp{—x/2 — 2)[1 + 0{l/x)] is small but greatly exceeds the value 
of E2 — El, e. g., Ei — E3 ^ 0.3 meV at i? = 60 nm. The donor asymmetry in the 
presence of the gate potentials will result in a further increase in E^ — E3, so that the 
phonon absorption processes will not contribute much to decoherence at sufficiently low 
temperatures T < 10 K. 

Thus, the error rate due to phonon- induced decoherence is Dijop) ~ 3 • 10~^ at Top < 
100 ps and T < 10 K. This value should be compared to the error rates due to other 
sources of decoherence. The lowest bounds for the decoherence times associated with the 
Johnson noise from the gates and the environmental charge fluctuations are fOJ 1201 120] , 
respectively, r ~ 1 /is and r ~ 1 ns, so that the corresponding error rates pHJ D{Top) = 
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1 — exp(— Top/r) do not exceed that due to phonons at Top < (1 ^ 10) ps. Hence, the 
performance of the buried donor charge qubit appears to be hmited primarily by the 
electron-phonon interaction. In this paper, we concentrated on the phosphorous donors 
in sihcon. Since the spectral function (jH^ that ultimately determines the error rate for 
one-qubit operations is a material constant, it would be worthwhile to search for other 
materials and/or doping elements for the buried donor charge qubit, in order to weaken 
the decoherence effects. 

Although we restricted ourselves to rectangular shapes of the resonant pulses, our 
consideration can be generalized to other pulse shapes [H]. The results obtained can be 
also applied to quantum-dot structures and Josephson three-level gates jlH ESI HI] • 
Finally, once a fundamental possibility of the auxiliary-state-assisted operations has been 
demonstrated, it is straightforward to organize the coupling of P2^:Si qubits for conditional 
quantum operations jTB} 1^. 

In summary, we have proposed a scheme for fast rotations of the buried donor charge 
qubit through an auxiliary-state-assisted electron evolution under the influence of resonant 
microwave pulses. This scheme allows for implementation of one-qubit operations in 
times as short as Top ~ 1 ps. By the example of the P^:Si qubit, we have shown that 
dephasing and "non- resonant" relaxation due to acoustic phonons are the main sources 
of decoherence. The error rate at T < 10 K and operation times Top = (1 ^ 10) ps is 
about 3 • 10~^, i. e., greater than the fault-tolerance threshold for quantum computation 
but sufficiently low to investigate the small-scale devices and thus to demonstrate the 
experimental feasibility of the scheme. 
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1. The logical states |0) = \L) and |1) = \R) of the buried donor charge qubit. 
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Fig. 2a. 

Fig. 2. One-electron wave functions of the two lowest states, IscXg (a) and 2pau (b), 
and the excited state 3dag (c) of the molecular ion P2':Si in the isotropic effective mass 
approximation. The coordinate X is along the line joining the two donors. The donor 
separation is i? = 20a^. The symmetric and antisymmetric linear superpositions of IsUg 
and 2p(Ju states correspond to Is atomic states \L) and \R) localized at the left and the 
right donor, respectively. They form the qubit logical states |0) = \L) and |1) = \R). 
The excited state 3dag is an auxiliary ("transport") state needed to transfer an electron 
between \L) and \R) states upon the influence of the external electromagnetic field. 
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Fig. 2b. 
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Fig. 2c. 
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